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1 Introduction 

Over the past few years the mathematical structure of the static correlation functions of the 
XXZ chain was largely resolved. After an appropriate regularization by a disorder parameter 
they all factorize into polynomials in only two functions p and uj [8]. These are the one-point 
function and a special neighbor two-point function which, in turn, can be represented as integrals 
over solutions of certain linear and non-linear integral equations [2]. This resembles much the 
situation with free fermions, and what is behind is indeed a remarkable fermionic structure on 
the space of quasi-local operators acting on the spin chain [5]. It allows us, for instance, to 
calculate short-range correlators with high numerical precision directly in the thermodynamic 
limit HKH]. 

The low temperature asymptotics of p and uj universally determines the low temperature 
properties of all static correlation functions. In this short note we obtain the low temperature 
asymptotics in the special case of the isotropic Hamiltonian 

with no magnetic field applied and vanishing disorder parameter. Then p = 1 and we are left 
with only one function (and its derivatives) which, up to a trivial factor, is the function 7 defined 
in [3]. 



2 Definition of the basic function 7 

For our purpose here it is convenient to introduce the function 7 within the context of a special 
realization of a six- vertex model (see e.g. [4]) and its associated quantum transfer matrix jlOj . 

*This paper is a contribution to the Proceedings of the International Workshop "Recent Advances in Quantum 
Integrable Systems" . The full collection is available at |http://www. emis.de/journals/SIGMA/RAQIS2010. html. 
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By definition the latter has 2{N vertical lines alternating in direction and carrying spectral 
parameters 

n, —u, u, —u, . . . , —u, u' + Hi, Hi — u', u' + fii, fj-i — u' , . . . , Hi ~ ■ 



2N 2M 

The spectral parameter on the horizontal line will be denoted H2- We consider this system in 
the limit J\f,A4 — )• +00 with the fine tuning uj\f = and u'Al = ly. With an appropriate 
overall normalization the largest eigenvalue A{h2,Hi) is given by 

AttJ 4:Tt5 
ln(A(/i2, ^1)) = -jrK{f^2) + —K{h2 - Hi) 

^ r ^^ In [{l + b{t,Himi + Ht,Hi))] ^2 1) 

y_oo 2cosh(7r(/i2 - i)) 

Let us note that we recover the familiar system of equations, allowing us to study the thermo- 
dynamical properties of the Hamiltonian by setting 6 = 0. The function K{x) is defined 

as 

1 ^~ikx 



1 

K(x) = — dk 



1 / , .x\ fl+ix\ fl-ix\ , ,x 

= TA*V-'2)-^\—)-*[—)^*i'^'-2 

where ■0 is the digamma function. The auxiliary functions b{x,H) a-nd h{x,H) ^I's solutions of 
a pair of non-linear integral equations given by 

ln(b(x, Hi)) = ^ - T ^ + r - + ^1)) 

1 cosh(7rxj 1 cos\i[tt[x — Hi)) J~oo 

/oo 
dtK{x-t + i)\n{l + b{t,Hi)) (2-2) 
-00 



and a similar equation obtained by exchanging b -f-)- b and i -f-)- — z in (|2.2p . The function 7 can 
now be introduced as 

7(^1,^2) = -! + (! + {Hi-H2f) r^ln(A(/i2,^i))|^^^. (2.3) 

It has been conjectured [3] that the correlation functions of the isotropic Heisenberg chain at 
any finite temperature (for vanishing magnetic field) are polynomials in 7 and its derivatives 
evaluated at (0,0). A similar statement (involving a function uj and its derivative with respect 
to the disorder parameter) was proved for the anisotropic XXZ chain O [HI [2] . Amazingly the 
isotropic limit seems non-trivial and is still a subject of ongoing work. Here we would only like 
to mention that the nearest- and next-to-nearest-neighbor two-point functions were expressed 
in terms of 7 in [3] starting from the multiple integral representation for the density matrix of 
the Heisenberg chain obtained in [7j. The formulae for the longitudinal two-point functions are, 
for instance, 

(cTfcT|)T = -^7(0,0), (2.4) 

{uIuDt = -^7(0,0) - ^7xx(0,0) + ^7xj;(0,0). (2.5) 

They will be used below to test our results for the low-temperature expansion. We denoted 
derivatives with respect to the first (resp. second) argument by the subscript x (resp. y). Similar 
results for four sites can be obtained from [1] in the isotropic limit. In previous work [13] the 
high-temperature expansion (up to order 25) of the two-point functions was obtained analytically 
based on ((231) and (1231). 
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3 Low-temperature expansion 

To compute the low-temperature expansion of 7, we follow the line of reasoning of the article [S], 
where a similar task was performed for the free energy. There are, however, two differences 
between the usual equations and the ones used in this note: the additional driving term in (j2.2p 
proportional to 5 and the shift /.f2 in the kernel of the integration in (j2.ip . 

The computation is based on the introduction of a shift C = -^In (vr^) in the auxiliary 
functions: 

hc{x) = b{x + C) and bc{x) = b{—x — C). 
In the low-temperature limit these functions satisfy 

Hbcix, ^1)) ~ -4e-"" - 4^ e-"("-'^i) + Vc{x) 

/oo 
dt[K{x - t) ln(l + bc{t, w)) -K{x-t + i) ln(l + bc{t, fii))] , (3.1) 

where Vc^x) is the rest of the integral which does not contribute to the low-temperature limit, 
when the magnetic field vanishes (see [9]). A similar relation holds with bob and i o — i 
exchanged. 

In terms of the shifted functions the largest eigenvalue becomes 

4tt J 4:Tr5 
ln(A(//2, /^i)) ~ -jrK{ji2) + -jrKin2 - f^i) 

+ dte''^^'^-'hn[{l + bcit,fii)){l + bcit,fii))] 

+ 7^ / dte-^^'^'+'hn[{l + bc{t,f,i)){l + bc{t,fii))]. 

J — c 

To evaluate these integrals we compute 

/oo 
^ dt[\ii{l + bc{t, W)) ln(b/:(t, + ln(l + bcit, /ii)) ln(b£(t, /ii))'] 

using two different methods. Here the prime stands for the derivative with respect to t. First, 
we compute it explicitly using the change of variables z = ln(b£) or z = ln(b£), respectively, 
which results in 



1 = 2 \n{l + e'')dz 
J —00 



7r2 
6 ■ 



Second, we replace ln(b£(t, //i)) and ln(b£(t, //i)) by their scaling limits (13. 1|) and simplify the 
resulting expression by taking into account that the derivative of K{x) is odd and contributions 
by double integrals cancel pairwise. This way we obtain 



X = 47r (^1 + -f""^^^ dte-""' In [(1 + b(t, + b(t, ^1))] . 



The same type of manipulation can be performed for the functions 6, and a similar result is 
obtained with fii replaced by —fii. 

Gathering these findings we obtain the asymptotic form of the largest eigenvalue, 

WA/ ^^ 47rJ^^, . 47r5^^, , T I e^^2 e-^^'^ \ 

ln(A(/i2, ^n)) ~ -7^Kifi2) + -7TrK{n2 - /^i) + tttt + 



24J I 1 + le'^A'i 1 + ie-'^^i 
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Figure 1. Comparison of the high- and low-temperature expansions (HTE, LTE) of {(jfa^) with the 
full numerical solution obtained from the integral equations (NLIE) and with Monte-Carlo data (QMC). 

Thus, using (|2.3p . the function 7 behaves asymptotically for small temperatures as 



This is our main result. 

Using (12. 4p and (12. Sp , we obtain the low-temperature expansion of the longitudinal correlation 
functions 



The constant terms (independent of the temperature) in these expansions are in agreement 
with those originally found in [11^ [6]. In the figure we compare the combined low- and high- 
temperature results for the next-to-nearest neighbor zz-correlation functions with the full nu- 
merical curve obtained by implementing the linear and non-linear integral equations that deter- 
mine 7 and its derivatives [3] on a computer. The high-temperature data and some additional 
Monte-Carlo data are taken from [14]. We find that the numerical curves (NLIE, QMC) are 
amazingly well approximated by its low- and high-temperature approximations. 
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7(a*i> /"2) ~ -1 + (1 (/Ui - fi2f) ( 4:TtK{h2 - m) 




W^I)T~^-^ln(2) + 
, ^ ^> 1 16 , , , 
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